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1. Introduction 

"In the thermodynamic limit the canonical and the microcanonical ensemble 
are equivalent in all details and generality. " Statements like this are found 
in many textbooks of statistical thermodynamics. It is the purpose of this 
contribution to show that this is not so, and, more importantly, that the 
microcanonical ensemble allows for significant insight into the mechanism 
of first order phase-transitions which is hidden in the canonical ensemble. 
E.g.: at the liquid-gas transition under given pressure large and fluctuating 
spatial inhomogeneities are created. This surface entropy S sur f governs the 
surface tension. At the critical point T cr S sur f compensates even the surface 
energy of the ground state leading to the vanishing of the surface tension. 

In the canonical ensemble at given pressure, however, inhomogeneities 
become suppressed by 

^ e -N 2 / 3 f surf /T bml ^ 

relative to the pure liquid or pure gas configurations, where N is the number 
of particles and f sur f the free energy of the interphase surface per surface 
atom.( In the following we use small letters for the energy or entropy per 
atom.) Consequently, macroscopic interfaces do not exist in the canonical 
ensemble. 

This has a very practical consequence: "If phase transitions in the exper- 
imental world would be at constant temperature not under controlled supply 
of energy" i.e. energy constraints unimportant as is often claimed, we would 
not be able to see a pot with boiling water. 
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Other common wisdoms are: "Phase transitions exist only in the ther- 
modynamic limit. " This is certainly true within the canonical ensemble 
where phase transitions are indicated by a singularity point as functions of 
the temperature. For finite systems, however, the canonical partition sum 



is analytic in (3 = 1/T. In the micro ensemble the transition occurs over a 
whole region in energy, relatively independent of the system size N. 

"S must be additive (extensive). For microscopic objects like atoms and 
cosmic like stars thermodynamics is inappropriate !"[l]In contrast, I will 
show in the following that phase transitions are well defined and can be 
well classified into first order or continuous transitions for pretty small 
systems like nuclei or some 100—1000 atoms in the microcanonical ensemble. 
Moreover, for the liquid-gas transition in sodium, potassium, and iron the 
microcanonical transition temperature Tj r , the specific latent heat qi a t, and 
the specific interphase surface tension a sur f are for ca. 1000 atoms close 
to their bulk values [2]. Nothing special qualifies the liquid-gas transition 
in infinite systems (at least for the above three realistic systems). This 
is also so in several toy-models of statistical mechanics: In [3] we found 
that for a microcanonical 2-dim Potts model with q = 10 states/lattice- 
point Tf r , qi a t , cr sur f are within a few percent close to their bulk values for 
some hundred spins. This astonishing fast scaling of the microcanonical 
transition parameters with the particle number towards their bulk values 
was explained by Hiiller and Promberger in a recent paper [4] by the fact 
that the trivial factor N in the exponent of the Laplace transform (2) is 
the main origin of canonical finite size scaling, which has nothing to do 
with the physics of the phase transition: Even if one replaces sjv(e) in (2) 
by Sqo (e) of the infinite system one gets practically the same scaling of the 
specific heat with N. 

In macroscopic thermodynamics extensive variables like total energy, 
mass, and charge are fixed only in the mean by the intensive T, fi^, \i% [5]. 
In nuclear physics they are, however, strictly conserved. Nuclear systems 
are too small to ignore the fluctuations of conserved quantities when using 
the usual macroscopic thermodynamic relations as Hill suggests even if the 
fluctuations are somehow ~ \ j\[N . Especially, if we are interested in phase 
transitions, we must be more careful: At phase transitions the effect of the 
fluctuations is usually even much larger, see below. Therefore, we test here 
statistical mechanics under extremely different conditions than we are used 
to. The foundations of thermodynamics must be revisited, extended and 
deepened. The constraints due to global conservation laws must carefully 
be respected and separated from the still persistent statistical fluctuations. 
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What is more remarkable and often overseen: Even in macroscopic sys- 
tems statistical fluctuations/parfide do not vanish and are usually large at 
phase-transitions of first order. I.e. even for phase transitions in the bulk 
(the standard field of thermodynamics) fluctuations and their constraints 
are important. Because it takes fluctuations serious our extended statisti- 
cal thermodynamics can be applied to phase transitions in small as well 
in large systems. Especially the development of phase transitions with the 
size can then be followed from small to the bulk. Microcanonical thermody- 
namics even sheds new light on details of the mechanism leading to phase 
transitions in the bulk (sect. 3). Moreover, in contrast to [1] it allows to ad- 
dress thermodynamically unstable like gravitating systems. There is some 
chance, this is possible even without invoking an artificial container. 

Another example illuminates what I want to emphasize: Angular mo- 
mentum is certainly one of the more exotic conservation laws in thermo- 
dynamics. In reality, however, it is much more important than one might 
think: In astrophysical systems like e.g. collapsing cosmic hydrogen clouds 
it may decide if a single star or a rotating multiple star system is born. 

2. Microcanonical thermodynamics 

Microcanonical thermodynamics explores the topology of the N-body phase 
space and determines how its volume W(E, N) = e s depends on the fun- 
damental globally conserved quantities namely total energy E = N * e, 
angular momentum L, mass (number of atoms N), charge Z, linear mo- 
mentum p, and last not least, if necessary, the available spatial volume V 
of the system. This definition is the basic starting point of any statisti- 
cal thermodynamics since Boltzmann[6]. If we do not know more about a 
complicated interacting N-body system but the values of its globally con- 
served macroscopic parameters, the probability to find it in a special phase 
space point (N-body quantum state) is uniform over the accessible phase 
space. This is an entirely mechanistic definition. It is of course a completely 
separated and difficult question, outside of thermodynamics, if and how a 
complicated interacting many-body system may explore its entire allowed 
phase space. In a nuclear or cluster collisions it is not really necessary that 
every dynamical path visits the entire possible phase space. It is sufficient 
that the evolution of an ensemble of many replica, one after the other, of 
the same system under identical macroscopic initial conditions follows the 
structure of the underlying N-body phase space. It is ergodic in the same 
sense as the dynamics of a falling ball is ergodic on Galton's nailboard. In 
nuclear fragmentation the ergodicity is presumably due to the strong and 
short ranged friction between moving nuclei in close proximity. Friction be- 
tween atomic clusters is yet unknown but probably it exists there also. 
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2.1) Statistical ensemble 

Before we proceed, we have to emphasize the concept of the statistical 
ensemble under strictly conserved energy, angular momentum, mass, and 
charge. As we explained above, microcanonical thermodynamics describes 
the dependence of the volume e s of the - at the given energy E, small in- 
terval SE, - accessible phase space on the globally conserved energy, mass, 
and charge. Each phase space cell of size (27rft) 3Af ~ 6 corresponds to an in- 
dividual configuration (event) of our system. We realize the ensemble by 
replica in time under identical macroscopic conditions (events) . This is dif- 
ferent from Hill [5] who assumes a macroscopic noninteracting supersystem 
of many identical copies of the system under consideration. This would be 
impossible e.g. for rotating or gravitating systems. Clearly, the volume e s 
of the phase space is the sum (ensemble) of all possible phase space cells 
compatible with the values of energy etc.. 

While the conserved, extensive quantities, energy, momentum, number 
of particles, and charge can be determined for each individual configura- 
tion of the system, i.e. at each phase-space point or each event, this is 
not possible for the phase space volume e s , i.e. the entropy S(E,V, N) 
and all its increments like the temperature T = (dS(E,V,N)/dE)~ 1 , 
the pressure P(E,V,N) = TdS(E,V, N)/dV , and the chemical potential 
fi = —TdS(E, V, N)/dN. They are ensemble averages. Only in the thermo- 
dynamic limit N — > oo may e.g. the temperature be determined in a single 
configuration by letting the energy flow into a small thermometer. For a 
finite system, e.g. a finite atomic cluster, the temperature, its entropy, its 
pressure can only be determined as ensemble averages over a large number 
of individual events. E.g. in a fusion of two nuclei the excitation energy in 
each event is given by the ground-state Q-values plus the incoming kinetic 
energy whereas the temperature of the fused compound nucleus is deter- 
mined by measuring the kinetic energy spectrum of decay products which 
is an average over many decays. It is immediately clear that the size of S 
is a measure of the fluctuations of the system. 

2.2) Microcanonical signal for a phase transition: The caloric curve. 
The most dramatic phenomena in thermodynamics are phase transitions. I 
will try to interpret them micro canonically as peculiarities of the topology 
of the N-body phase space. I will avoid the concept of the thermodynamic 
limit as I believe that this is not really essential to understand phase tran- 
sitions. We will see that details about the transitions become more trans- 
parent in finite systems. Then however, one needs a modified definition of 
phase transitions. 

In [7, 8] we introduced a new criterion of phase transitions, which 
avoids any reference to the thermodynamic limit and can also be used 
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for finite systems: The anomaly of the micro canonical caloric equation of 
state T(E/N) where dT/dE < i.e. where the familiar monotone rise 
of the temperature with energy is interrupted. Rising the energy leads 
here to a cooling of the system. This anomaly corresponds to a convex 
intruder in S(E) = J l/T(E)dE. At energies where S(E) is convex the 
system would spontaneously divide into two parts and gain entropy eg.: 
+ S{E 2 ))/2 > S((E X + E 2 )/2), iff the creation of the interface 
would not cost an extra entropy AS sur f. Because As sur f per atom van- 
ishes oc N^ 1 / 3 in the thermodynamic limit Soo(e) is concave as demanded 
by van Hove's theorem [9] and by the second law of thermodynamics. 

Very early the anomaly of the caloric curve T(E/N) was taken as signal 
for a phase transition in small systems in the statistical theory of multi- 
fragmentation of hot nuclei by Gross and collaborators [7, 10] and the 
review article [11]. Bixton and Jortner [12] linked the back-bending of the 
microcanonical caloric curve to strong bunching in the quantum level struc- 
ture of the many-body system i.e. a rapid and sudden opening of new phase 
space when the energy rises. Their paper offers an interesting analytical in- 
vestigation of this connection. 

A phase transition of first order is characterized by a sine-like oscil- 
lation, a "back-bending" of T(e = E/N) c.f. fig.l. As shown below, the 
Maxwell-line which divides the oscillation of dS/dE = (3(e) = 1/T into 
two opposite areas of equal size gives the inverse of the transition temper- 
ature Tf r , its length the specific latent heat qi a t, and the shaded area under 
each of the oscillations is the loss of specific entropy As sur f as mentioned 
above. The latter is connected to the creation of macroscopic interphase 
surfaces, which divide mixed configurations into separated pieces of dif- 
ferent phases, e.g. liquid droplets in the gas or gas bubbles in the liquid. 
Even nested situations are found like liquid droplets inside of crystallized 
pieces which themselves are swimming in the liquid in the case of the solid 
- liquid transition, see e.g. the experiments reported in [13]. I.e. at phase 
transitions of first order inhomogeneous "macroscopic or collective" density 
fluctuations are common, boiling water is certainly the best known exam- 
ple. Phase-dividing surfaces of macroscopic size exist where many atoms 
collectively constitute a boundary between two phases which cause the re- 
duction of entropy by As sur f. 

2.3) "Maxwell" construction ofTt r , qi a t, o~ sur f 

As the entropy is the integral of (3(e): s(e) = Jq (3(e')de' it is a concave 
function of e (d 2 s/de 2 = d(3/de < 0) as long as T(e) = (3" 1 shows the 
usual monotonic rise with energy. In the pathological back-bending region 
of (3(e) the entropy s(e) has a convex intruder of depth As sur j [3]. At the 
beginning (> ei) (c.f. fig.lc) of the intruder the specific entropy s(e) is re- 



6 



D.H.E. GROSS 



duced compared to its concave hull, which is the double tangent to s(e) in 
the points e\ and e%. The derivative of the hull to s(e) follows the Maxwell- 
line in the interval e\ < e < e^. In the middle, (e2), when the separation 
of the phases is fully established this reduction is maximal = As sur f and 
at the end of the transition (e^) when the intra-phase surface(s) disappears 
As sur f is gained back. Consequently, the two equal areas in /3(e) are the 
initial loss of surface entropy As sur f and the later regain of it. Due to van 
Hove's theorem this convex intruder of s(e) must disappear in the thermo- 
dynamic limit which it will do if As sur f ~ N^ 1 / 3 . This is why a transition 
of first order may easier be identified in finite systems where the intruder 
can still be seen. The intra-phase surface tension is related to As sur f by 
7 sur/ = As SMr/ * N * T tr /surf. -area. [14]}. 




0.3 0.5 0.7 0.9 1.1 1.3 



Excitation Energy (eV/Atorn) 
Fig.l: microcanonical caloric curve T P (E/N = e,V(E, P, N), N) = (dS/dE)~ 1 \ P at constant 
pressure (full square points), number of fragments Nf r with m; > 2 atoms, and the effective 
number of surface atoms JV^yy = ^2 n h^' = total surface area divided by 4irr^ ua . In the lower 
panel s(e) = f e f3(e')de' is shown. In order to make the intruder between e\ and eg visible, we 
subtracted the linear function 25 + 11. 5e. 

2.4) Geometric vs. entropic classification phases and phase-boundaries 
Our approach to phase-transitions of first order is complementary to the 
conventional approach where the separation of the system into two homo- 
geneous phases by a — in general — geometrical interface is investigated, 
e.g. [15]. The problems due to the large fluctuations of this interface are 
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numerous and severe, see e.g. the discussion on fluid interfaces by Evans 
[16]. These fluctuation are of course crucial for the interfacial entropy and 
consequently for the surface tension also. 

The two main differences of our approach are that a) we do not start 
with the geometry (planar or spherical) of the interface, nor do we demand 
the interface to compact. We focus our attention to the entropy of the 
phase separation. This turns out to be much simpler than the geometric 
approach, b) The microcanonical ensemble allows for large scale spatial 
inhomogeneities, whereas the canonical ensemble suppresses spatial inho- 
mogeneities like phase-separations exponentially. Any interface needs an 
additional free energy. {In the case of a phase transition of first order the 
suppression is oc exp(—aN 2 / 3 /T tr ), where a is the surface tension parame- 
ter (a = 47rr^ s 7) and r ws is the radius of the Wigner-Seitz cell. 

Our characterization of phase transitions is purely thermodynamically. 
We have not yet defined what a phase is. Much effort is spent by Ruelle to 
define pure phases as those for which in the thermodynamic limit observ- 
ables survive increasing coarse-graining, for which space- averaged quanti- 
ties do not fluctuate, c.f. chapter 6.5 in Ruelle's book [17]. Of course, this 
definition works in the thermodynamic limit only. It does not address finite 
systems. For a finite system it is not possible to decide if a single config- 
uration corresponds to a pure phase or not. The situation is analogous to 
the definition of the temperature, see above. Again to be a pure phase is 
a feature of the whole ensemble not of a single phase-space point (config- 
uration). We offer a statistical definition of a pure phase: A configuration 
belongs to the ensemble of pure phases — including its fluctuations — at 
concave points of S(E) with d 2 S/dE 2 < 0. 

2.5) The microcanonical constant pressure ensemble 
The micro-canonical ensemble with given pressure {E,P,N} must be dis- 
tinguished from the (in spirit) similar constant pressure ensemble {H,P,N} 
introduced by Andersen [18, 19, 20] where a molecular-dynamic calculation 
with the hypothetical Hamiltonian [19] 

T/ 2 / 3 N N N M 

H({ ri , Pi }, V,V) = -—^2 mn 2 + £ ]T Hr^V 1 / 3 ) + —V 2 + P E V (3) 

i=i j=ij>i 

is suggested. Here V is the volume of the system, taken as an additional 
explicit degree of freedom, {ri, P i} are the coordinates and momenta of the 
atoms scaled with the factor F 1 / 3 , Ti,V the corresp. velocities, $(r^) is 
the intra-atomic two-body potential, and M is a hypothetical mass for the 
volume degree of freedom. Pe is the given pressure. The total enthalpy H, 
atoms plus V-degree of freedom, is conserved, not the total energy E of the 
atoms alone. 



8 



D.H.E. GROSS 



This is very different to our micro-canonical approach with given E, 
V(E, P) ,N where the energy E of the atoms is conserved and the pressure 
is the correct thermodynamic pressure (P(E,V) = T(E,V)dS/dV\E)- At 
each energy the volume V(E, P) is simultaneously chosen for all members 
of the ensemble by the condition that T(E, V)dS(E, V)/dV\ E of the whole 
ensemble is the correct pressure. In this case there is a unique correlation 
between the energy E and the volume V which does not fluctuate within 
the ensemble even though the pressure is specified. At the given energy this 
is still the {E,V(E,P,N),N} ensemble. 



3. The liquid-gas transition of sodium, potassium, and iron 

The microscopic simulation of the liquid-gas transition in metals is espe- 
cially difficult. Due to the delocalisation of the conductance electrons metals 
are not bound by two-body forces but experience long-range many-body in- 
teractions. Moreover, at the liquid-gas transition the binding changes from 
metallic to covalent binding. This fact is a main obstacle for the conven- 
tional treatment by molecular dynamics [21]. 

In the macro-micro approach we do not follow each atom like in molec- 
ular dynamics, the basic particles are the fragments. Their ground-state 
binding energies are taken from experiments. The fragments are spherical 
and have translational, rotational, and intrinsic degrees of freedom. The 
internal degrees of freedom of the fragments are simulated as pieces of bulk 
matter. The internal density of states, resp. the internal entropy of the 
fragments is taken as the specific bulk entropy s(e) at excitation energies 
e < e max = e-boih which can be determined from the experimentally known 
specific heat of the solid/liquid bulk matter [22]. e^n is the specific energy 
where the boiling of bulk matter starts. Details are discussed in [2, 14]. Then 
the metallic binding poses no difficulty for us and the metal — nonmetal 
transitions is in our approach controlled by the increasing fragmentation 
of the system. This leads to a decreasing mean coordination number when 
the transition is approached from the liquid side while the distance to the 
nearest neighbor keeps about the same. Exactly this behavior was recently 
experimentally observed [23, 24]. By using the microcanonical ensemble we 
do not prespecify the intra-phase surface and allow it to take any form. 
Also any fragmentation of the interface is allowed. It is the entropy alone 
which determines the fluctuations of the interface. Here we present the first 
microscopic calculation of the surface tension in liquid sodium, potassium, 
and iron. 

The decay of potassium is in all details similar to that of sodium, fig.(l). 
Therefore we don't show here the corresponding figures. The liquid-gas 



PHASE TRANSITIONS WITHOUT THERMODYNAMIC LIMIT 



9 



transition in iron is different from that of the alkali metals: Due to the 
considerably larger surface energy parameter a s in the liquid drop formula 
of the ground-state binding energies of iron compared to alkali metals there 
is no multi-fragmentation of iron clusters at P = 1 atm. Iron clusters of 
N < 3000 atoms decay by multiple monomer evaporation c.f. fig.2. 
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TABLE llParameters of the liquid— gas transition at constant pressure of 1 atm. in a mi- 
crocanonical system of No interacting atoms and in the bulk. st, n = qi a t/Ttr, it is interesting 
that the value of s^ a for all three systems and at all sizes is near to s^oU = 10 as proposed 

by the empirical Trouton's rule [25], N^tt = m^/ 3 Nj, Nf = multiplicity of fragments with 

2/3 

mi-atoms, and cr/Ttr = NoAs sur f /N Jff The bulk values a/Ttr are adjusted to the input values 
of a s taken for the T = surface tension from ref.[26] which we used in the present calculation 
for the ground-state binding energies of the fragments. Ttr is Tp not Ty. Ty comes even closer 
to the bulk value. Other inputs of the calculations are: gound state liquid drop parameters of all 
possible clusters and the bulk entropy/internal dof in the liquid phase at e < ei (fig.l). 
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Table (1) gives a summary of all theoretical parameters for the liquid- 
gas transition in clusters of Nq = 200 — 3000 Na, K, and Fe atoms and 
compared with their experimental bulk values. The transition-temperature 
T tr , the specific latent heat qi a t and the entropy gain of an evaporated atom 
Sbou are approaching the experimental bulk values. As sur f is the area under 
the back-bend of /3(e). NoAs sur f is the total entropy loss due to the inter- 

^ 2/3 2/3 

faces equal to J2^ r ws m i ^ill^tr = ^ e ff a /Ttr- Of course, the transition 
temperature T tr and the latent heat qi a t of small clusters are smaller than 
the bulk values because the average coordination number of an atom at the 
surface of a small cluster is smaller than at a planar surface of the bulk. 



4. Conclusion 

We presented a formulation of statistical thermodynamics entirely from the 
principles of mechanics. By carefully distinguishing conserved quantities 
like the energy, or angular momentun from ensemble averaged quantities 
like entropy, temperature or pressure this formulation allows an application 
to small systems like nuclei or atomic clusters and to astrophysical systems. 

The microcanonical statistics also gives a detailed insight into the na- 
ture of phase transitions. The liquid-gas transition in metals at normal 
pressure of 1 atm. is experimentally well explored. Therefore, it is a good 
test case for our ideas and computational methods of microcanonical ther- 
modynamics. By allowing a system of N = 200 — 3000 atoms to condense 
or fragment into an arbitrary number of spherical fragment clusters and 
into an arbitrary number of free atoms under a prescribed external pres- 
sure of 1 atm. we calculated by microcanonical Monte Carlo methods [14] 
the microcanonical caloric curve (3{e) = ds/de =< d/de >. The anomaly 
of T(e), where dT/de < 0, signals the liquid-gas phase-transition in the 
finite system. The characteristic parameters T tr , qiat-, and the surface ten- 
sion a sur f are for ~ 1000 atoms smaller but similar to the experimentally 
known values of the bulk liquid-gas transition. 

This result is remarkable for several reasons: 

a) It proves that a phase transition of first order in a realistic continuous 
system can very well be classified in small mesoscopic clusters without in- 
voking the thermodynamic limit. In fact, with the finite back-bending of 
T{E/N) the transition is easier recognizable than at N — > oo. 

b) Consistent with similar results for statistical toy models [3] T tr , qiat, o~ surf 
reflect for astonishing small systems their bulk values. The mechanism lead- 
ing to phase-transition has nothing to do with the thermodynamic limit. 

c) Even for a realistic metallic system with its long-range many-body in- 
teractions the Mcrocanoncal Metropolis Monte Carlo simulation method is 
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able to describe the liquid-gas phase transitions quite well. This is possible 
because we do not use molecular dynamics with a two-body Hamiltonian 
but use the experimental groundstate binding energies of the fragment clus- 
ters which of course take care of the metallic bonding of their constituents. 

d) The intra-phase surface entropy can be microscopically calculated. This 
was not possible up to now for metals. The surface tension per surface 
area can be determined if the intra-phase area is known. For the surface 
entropy the fluctuation and fragmentation of the surface are essential. At 
the liquid-gas transition of sodium and potassium clusters of sizes as con- 
sidered here at normal pressure the intra-phase fluctuations are mainly due 
to strong inhomogeneities and clusterization. This is consistent with recent 
experimental evidence for the bulk [23, 24]. This result is encouraging to 
also compute the surface tension at higher pressure to see its vanishing to- 
wards the critical point. 

e) It is promising that Microcanoncal Metropolis Monte Carlo obtains val- 
ues for small clusters of comparable size to the known infinite matter values 
of the liquid-gas transition. This is a necessary test of MMMC to describe 
nuclear fragmentation [11, 14] and to be able to get insight into the (criti- 
cal) behavior of the nuclear matter liquid-gas transition. This is important 
as it is obviously not possible to compare the model with experimental data 
for nuclear matter. 

f) The division of the macroscopic observables of the N-body system into 
observables which can be determined at each phase-space point, in each in- 
dividual realization of the microcanonical ensemble, i.e.: globally conserved 
quantities, energy, number of particles, charge, angular momentum, 

and into thermodynamic observables which refer to the size of the ensem- 
ble, the volume e s of the energy shell of the phase space, and which cannot 
be determined at a single phase-space point, in a single event, like entropy, 
temperature, pressure, chemical potential is very essential. The concept of 
a pure phase and of a phase-transition belongs to the second group. 

g) Last not least, there is a very fundamental difference in the microcanon- 
ical and the canonical ensembles. The energy is the crucial parameter that 
controls the fluctuations in their development from the liquid to the vapor 
and this is respected only in the micro-ensemble. 

I am grateful to Th. Klotz, V. Laliena, O. Fliegans, M. Madjet, and E. 
Votyakov for many numerical computations and helpful discussions. 
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